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We study the quantum kicked rotator in the classically fully chaotic regime K — 10 and for various 
values of the quantum parameter k using Izrailev’s A'^-dimensional model for various N < 3000, which 
in the limit N ^ oo tends to the exact quantized kicked rotator. By numerically calculating the 
eigenfunctions in the basis of the angular momentum we find that the localization length £. for fixed 
parameter values has a certain distribution, in fact its inverse is Gaussian distributed, in analogy 
and in connection with the distribution of finite time Lyapunov exponents of Hamilton systems. 

However, unlike the case of the finite time Lyapunov exponents, this distribution is found to be 
independent of N, and thus survives the limit N = oo. This is different from the tight-binding 
model of Anderson localization. The reason is that the finite bandwidth approximation of the 
underlying Hamilton dynamical system in the Shepelyansky picture (D.L. Shepelyansky, Phys. Rev. 

Lett. 56, 677 (1986)) does not apply rigorously. This observation explains the strong fluctuations 
in the scaling laws of the kicked rotator, such as e.g. the entropy localization measure as a function 
of the scaling parameter A = C/N, where £ is the theoretical value of the localization length in 
the semiclassical approximation. These results call for a more refined theory of the localization 
length in the quantum kicked rotator and in similar Floquet systems, where we must predict not 
only the mean value of the inverse of the localization length £ but also its (Gaussian) distribution, 
in particular the variance. In order to complete our studies we numerically analyze the related 
behavior of finite time Lyapunov exponents in the standard map and of the 2x2 transfer matrix 
formalism. This paper is extending our recent work (T. Manos and M. Robnik, Phys. Rev. E 87, 

062905 (2013)). 
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I. INTRODUCTION 

Time-periodic (Floquet) quantum systems, whose clas¬ 
sical analog is fully chaotic and diffusive, typically exhibit 
dynamical localization [Il,i , if a certain semiclassical 
condition is satisfied, as explained below. We study the 
periodically kicked rotator in the classically fully chaotic 
regime K = 10 using Izrailev’s fV-dimensional model 

for various N < 3000, which in the limit N ^ oo tends 
to the quantized kicked rotator. We restrict our analy¬ 
sis to the case K = 10 because this is empirically the 
most typical uniformly chaotic regime, apparently free of 
any islands of stability or acceleration modes Q- Due 
to the finiteness of N the observed (dimensionless) local¬ 
ization length of the eigenfunctions in the space of the 
angular momentum quantum number does not possess a 
sharply defined value, but has a certain distribution in¬ 
stead. Its reciprocal value is almost Gaussian distributed. 
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This might be expected on the analogy with the finite 
time Lyapunov exponents in the Hamiltonian dynamical 
systems. In order to corroborate the theoretical findings 
on this topics we perform in Secs. m and El the nu¬ 
merical analysis of the finite time Lyapunov exponents 
in the standard map (classical kicked rotator), especially 
the decOT of the variance. Indeed, in the Shepelyansky 
picture [8| the localization length can be obtained as the 
inverse of the smallest positive Lyapunov exponent of a 
finite 2fc-dimensional Hamilton system associated with 
the band matrix representation of the quantum kicked 
rotator, where k is the quantum kick parameter (to be 
precisely defined below). In this picture, N plays the 
role of time. However, unlike the chaotic classical maps 
or products of transfer matrices in the Anderson tight- 
binding approximation, where the mean value of the fi¬ 
nite time Lyapunov exponents is usually equal to their 
asymptotical value of infinite time and the variance de¬ 
creases inversely with time, as we also carefully checked 
(see Secs. El and IVI|) , here the distribution is found to be 
independent of It has a nonzero variance even in the 
limit N ^ oo. The reason is that the quantum kicked 
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rotator at = oo cannot be exactly modeled with finite 
bandwidth (equal to 2k) band matrices, but only approx¬ 
imately, such that the underlying Hamilton system of 
the Shepelyansky picture has a growing dimension with 
TV, implying asymptotically an infinite set of Lyapunov 
exponents and behavior different from the finite dimen¬ 
sional Hamiltonian systems. The observation of the dis¬ 
tribution of the localization length around its mean value 
with finite variance also explains the strong fluctuations 
in the scaling laws of the kicked rotator, such as e.g. the 
entropy localization measure as a function of the theoret¬ 
ical scaling parameter A, to be discussed below. On the 
other hand, the two different empirical localization mea¬ 
sures, namely the mean localization length as extracted 
directly from the exponentially localized eigenfunctions 
and the measure based on the information entropy of 
the eigenstates, are perfectly well linearly connected and 
thus equivalent. Therefore these results call for a refined 
theory of the localization length in the quantum kicked 
rotator and similar systems, where we must predict not 
only the mean value of the inverse localization length but 
also its (Gaussian) distribution, in particular the vari¬ 
ance. This paper is a follow-up paper of our recent work 
0 (Manos and Robnik 2013). 

The time-independent and time-periodic systems have 
much in common when discussing the localization prop¬ 
erties of the chaotic eigenstates and of the corresponding 
energy spectra. The main result of stationary quantum 
chaos (or wave chaos) [IIJ,® is the discovery that in 
classically fully chaotic, ergodic, autonomous Hamilton 
systems with the purely discrete spectrum the fluctua¬ 
tions of the energy spectrum around its mean behavior 
obey the statistical laws described by the Gaussian Ran¬ 
dom Matrix Theory (RMT) [l^, [l^, provided that we 
are in the sufficiently deep semiclassical limit. The latter 
semiclassical condition means that all relevant classical 
transport times are smaller than the so-called Heisen¬ 
berg time, or break time, given by tn = 2nh/ AE, where 
h = 27rfi is the Planck constant and AE is the mean 
energy level spacing, such that the mean energy level 
density is p{E) = 1/AE. This statement is known as 
the Bohigas - Giannoni - Schmit (BGS) conjecture and 
goes back to their pioneering paper in 1984 llJ . although 
some preliminary ideas were published in |14| . Since 
AE oc ?i-^, where / is the number of degrees of freedom 
(= the dimension of the configuration space), we see that 
for sufficiently small h the stated condition will always be 
satisfied. Alternatively, fixing the h, we can go to high 
energies such that the classical transport times become 
smaller than tu- The role of the antiunitary symme¬ 
tries that classify the statistics in terms of GOE, GUE or 
GSE (ensembles of RMT) has been elucidated in [T^, see 
also [i^ , and [l|, 0, [13, [ill • The theoretical foundation 
for the BGS conjecture has been initiated first by Berry 

a , and later further developed by Richter and Sieber 
, arriving finally in the almost-final proof proposed 
by the group of F. Haake (l9[ - [^ . 

Here it must be emphasized again that considering 


the chaotic eigenstates and their dynamical localiza¬ 
tion properties there are strong analogies between the 
time-periodic systems (like the kicked rotator) and time- 
independent systems (like static billiards) [l^l , where the 
Brody distribution [^ plays a key role, as discussed 
in 1,12311 . 

The paper is organized as follows: In Sec. [Til we de¬ 
fine the model, in Sec. m we present the evidence for 
and the description of the distribution of the localization 
measures, in Sec. HYi we study the finite time Lyapunov 
exponents of the classical standard mapping as a generic 
example of a chaotic area preserving mapping, in Sec. [Vj 
we study the finite time Lyapunov exponents of the prod¬ 
uct of two-dimensional random symplectic matrices de¬ 
scribing the tight-binding model of Anderson localiza¬ 
tion, in Sec. EH we present the high precision numerical 
results about the decay of the variance of the distribu¬ 
tion of the finite Lyapunov exponents of Sec. llVIandIVl 
and in Sec. IVHI we conclude and discuss the results in 
the broader theoretical perspective. 


II. THE KICKED ROTATOR, THE IZRAILEV 
MODEL AND THE DYNAMICAL 
LOCALIZATION 

The kicked rotator was introduced by Casati, Chirikov, 
Ford and Izrailev in 1979 [sO]. Here we follow our nota¬ 
tion Q. The Hamiltonian function is 

H =^ + VodTit) cos9. (1) 

Here p is the (angular) momentum, I the moment of 
inertia, Vq is the strength of the periodic kicking, 9 G 
[0, 27r) is the (canonically conjugate, rotation) angle, and 
Srit) is the periodic Dirac delta function with period 
T. Between the kicks the rotation is free, and thus the 
dynamics can be reduced to the standard mapping, 

T 

Pn+l — Pn E Vq Sin^yi_|_i, ^n-t-1 — 9ji -f —Pnj (2) 

as introduced in [3ll - [^ . The quantities {9n,Pn) refer 
to their values just immediately after the n-th kick. By 
using new dimensionless momentum //, we get 

Pn+l = Pn + K sin0„_|_i, 0„+l = 9n + Pn, (3) 

where the system has now a single classical dimensionless 
control parameter K = VqT jl. 

The quantum kicked rotator (QKR) is the quantized 
version of Eq. o, namely 

^ =+ cos6l. (4) 

The Eloquet operator F acting on the wavefunctions 
(probability amplitudes) ^>(9), 9 G [0, 27r), upon each pe- 
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riod (of length T) can be written as (see e.g. [I|, Chapter 

4) 

f iVo f ihT d^\ 

where now we have two dimensionless quantum control 
parameters 


k = 


h ’ 



( 6 ) 


which satisfy the relationship K = kr = VqT/I, K being 
the classical dimensionless control parameter of Eq. ©. 
By using the angular momentum eigenfunctions 

{0\n) = a„(6») = ^=exp(zn0), (7) 

V 27r 


where n is any integer, we find the matrix elements of E, 
namely 


Emn = (m|F|n) = exp ’"J„_m(fc), (8) 

where J^{k) is the i/-th order Bessel function. For a 
wavefunction 'ip{9) we shall denote its angular momen¬ 
tum component (Fourier component) by 


Un 


(<2tt 


{n\il))= / a*„{e)ijj{e) dO = 


1 


V^. 


ip{6) exp{—in9) d9. 


(9) 


The QKR has very complex dynamics and spectral prop¬ 
erties. As the phase space is infinite (cylinder), p S 
(—oo,-foo), 0 S [0,27r), the spectrum of the eigenphases 
of F, denoted by (j)n, or the associated quasienergies 
hun = fi(j)n/T, introduced by Zeldovich [s^, can be con¬ 
tinuous, or discrete 

The asymptotic localized eigenstates are exponentially 
localized. The (dimensionless) theoretical localization 
length in the space of the angular momentum quantum 
numbers is given below, and is equal (after introducing 
some numerical correction factor a^) to the dimension¬ 
less localization time tioc [Em (IT^ .given below]. We 
denote it unlike in reference @ and Q by £. Therefore, 
an exponentially localized eigenfunction centered at m in 
the angular momentum space [Eq. (O] has the following 
form 


^exp 



( 10 ) 


where is the probability amplitude [Eq. (0] of the lo¬ 
calized wavefunction ip{9). The argument leading to tioc 
in Eq. m given below originates from the observation 
of the dynamical localization by Casati et al 30| , and in 
particular from [s^, and is well explained in [ij, in case of 


normal diffusion, whilst for general anomalous diffusion 
we gave a theoretical argument in We shall denote 
a = 2/£, and will later on determine the cr’s directly 
from the individual numerically calculated eigenstate. 

The question arises, where do we see the phenom¬ 
ena (spectral statistics, namely Brody-like level spacing 
distribution) analogous in the quantum chaos of time- 
independent bound systems with discrete spectrum? To 
see these effects the system must have effectively finite 
dimension, because in the infinite dimensional case we 
simply observe Poissonian statistics. Truncation of the 
infinite matrix Fmn in Eq. ([8]) in tour de force is not 
acceptable, even in the technical case of numerical com¬ 
putations, since after truncation the Floquet operator is 
no longer unitary. 

The only way to obtain a quantum system which shall 
in this sense correspond to the classical dynamical sys¬ 
tem [Eqs. (fTI) . (1^ and ©] is to introduce a finite N- 
dimensional matrix, which is symmetric unitary, and 
which in the limit N ^ oo becomes the infinite dimen¬ 
sional system with the Floquet operator [Eq. ([5])]. The 
semiclassical limit is fc —>■ oo and r —> 0, such that 
K = kr = constant. As it is well known for the 
reasons discussed above, the system behaves very simi¬ 
larly for rational and irrational values of r/(47r). Such a 
Wdimensional model [1^ will be introduced below. 

The generalized diffusion process of the standard map 
© is defined by 

{iAPr)=D,{K)n>^, (11) 

where n is the number of iterations (kicks), and the ex¬ 
ponent p, is in the interval [0,2), and all variables P, 9 
and K are dimensionless. Here Df^{K) is the general¬ 
ized classical diffusion constant. The averaging (.) 
is over an ensemble of initial conditions with fixed P, 
specifically in our case P = 0. In case p = 1 we have 
the normal diffusion, and Di{K) is then the normal dif¬ 
fusion constant, whilst in case of anomalous diffusion we 
observe subdiffusion when 0 < /r < 1, or superdiffusion 
a 1 < p < 2. In case p = 2 we have the ballistic trans¬ 
port which is associated with the presence of accelerator 
modes (see below). 

Following 0 we find that the dimensionless Heisen¬ 
berg time, also called break time or localization time, 
denoted by fioc, in units of kicking period T, is equal to 
the dimensionless localization length C 


C 



1 

2 — p, 


( 12 ) 


where is a numerical constant to be determined em¬ 
pirically, and in case of normal diffusion p = I is close to 
1 / 2 . 

In case of the normal diffusion p = 1, considered in the 
present paper, the theoretical value of Di{K) is given in 
the literature, e.g. in @ or (4ll| . 
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D U^M1-2J2W(1-J2W)],ifi^>4.5 

’ [0.15(i^:-^Ccr)^if i^cr- < < 4.5 

(13) 

where Kcr — 0.9716 and J 2 {K) is the Bessel func¬ 
tion. Here we neglect higher terms of order K~^. In 
the present paper we shall consider exclusively the case 
K = 10, which has been carefully checked to be fully 
chaotic, without any regular islands, and well described 
by the normal diffusion /r = 1, so that the above formula 
applies very well Q- 

The motion of the QKR [Eq. (j?])] after one period T of 
the Ip wavefunction can be described also by the following 
symmetrized Floquet mapping, describing the evolution 
of the kicked rotator from the middle of a free rotation 
over a kick to the middle of the next free rotation, as 
follows 

pj{e,t + T) = u^p{e,t), (14) 

TT (Thd^\ ( .Vo A (.Thd^\ 

Thus, the ipiO, t) function is determined in the middle of 
the rotation, between two successive kicks. The evolution 
operator U of the system corresponds to one period. 

In the case K = kr ^ 1 the motion is well known to 
be strongly chaotic, for K = 10 certainly without any 
regular islands of stability, and also there are no acceler¬ 
ator modes, so that the diffusion is normal (^ = 1). We 
have carefully checked that the case X = 10 is the clos¬ 
est to the normal diffusion /r = 1 for all K G [0, 70]. The 
transition to classical mechanics is described by the limit 
/c —>■ 00 , r —7> 0 while K = const. We shall consider the 
regimes on the interval 3 < fc < 20, but will concentrate 
mostly on the semiclassical regime k > K, where r < 1. 

In order to study how the localization affects the sta¬ 
tistical properties of the quasienergy spectra, we use the 
model’s representation in the momentum space with a 
finite number N of levels [l-@, , which we refer to as 

Izrailev model 


N 

Un{t + T) = ^ UnmUmit), 71, m = 1, 2, ..., N . (15) 

m—1 

The finite symmetric unitary matrix Unm determines the 
evolution of an fV-dimensional vector, namely the Fourier 
transform u„(t) of "0(0, t), and is composed in the follow¬ 
ing way 

Unm — ^ ) Gnm'Un'm'Gn'mj (16) 

n'm' 

where Gw = exp (irU /4^ Sw is a diagonal matrix corre¬ 
sponding to free rotation during a half period T/2, and 


the matrix Bn'm' describing the one kick has the follow¬ 
ing form 


Br, 


1 


2N + 1 


2N+1 

E 

1=1 


X exp 


cos 


{n' — m') 


2t:1 


—ik cos 


2N+1 
2ttI 


— cos 


{n' + m') 


2ttI 


2N + I 


2N + I 


(17) 


The Izrailev model in Eqs. (HMni) with a finite number 
of states is considered as the quantum analogue of the 
classical standard mapping on the torus with closed mo¬ 
mentum p and phase 6, where Unm describes only the 
odd states of the systems, i.e. 'tp{9) = —'ip{—9), pro¬ 
vided we have the case of the quantum resonance, namely 
T = Attt/{2N + 1), where r is a positive integer. The ma¬ 
trix (1171) is obtained by starting the derivation from the 
odd-parity basis of sin(n0) rather than the general angu¬ 
lar momentum basis exp(m0). 

Nevertheless, we shall use this model for any value of 
T and fc, as a model which in the resonant and in the 
generic case (irrational T/(47r)) corresponds to the clas¬ 
sical kicked rotator, and in the limit N ^ oo approaches 
the infinite dimensional model [Eq. (I14p ]. restricted to 
the symmetry class of the odd eigenfunctions. It is of 
course just one of the possible discrete approximations 
to the continuous infinite dimensional model. 

The difference of behavior between the generic case 
and the quantum resonance shows up only at very large 
times, which grow fast with {2N -1-1), as explained in 
0. It turns out that also the eigenfunctions and the 
spectra of the eigenphases at finite dimension N of the 
matrices that we consider do not show any significant 
differences in structural behavior for the rational or ir¬ 
rational r/(47r), which we have carefully checked. In¬ 
deed, although the eigenfunctions and the spectrum of 
the eigenphases exhibit sensitive dependence on the pa¬ 
rameters T and fc, their statistical properties are stable 
against the small changes of r and fc. This is an advan¬ 
tage, as instead of using very large single matrices for 
the statistical analysis, we can take a large ensemble of 
smaller matrices for values of r and fc around some cen¬ 
tral value of T = To and fc = fco, which greatly facilitates 
the numerical calculations and improves the statistical 
significance of our empirical results. Therefore our ap¬ 
proach is physically meanin gfu l. Similar approach was 
undertaken by Izrailev (see [6| and references therein). 
In Fig. 1 of paper Q we show the examples of strongly 
exponentially localized eigenstates by plotting the natu¬ 
ral logarithm of the probabilities Wn = junp versus the 
momentum quantum number n, for two different matrix 
dimensions N. By calculating the localization length £ 
from the slopes tr = 2/C of these eigenfunctions using 
Eq. (fTOl) we can get the first quantitative empirical local¬ 
ization measure to be discussed and used later on. The 
new finding of this paper is that a has a distribution, 
which is close to the Gaussian (but cannot be exactly 
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that, because cr is a positive definite quantity). It does 
not depend on N and survives the limit N ^ oo. There¬ 
fore also C has a distribution whose variance does not 
vanish in the limit iV —^ oo. 

Following 0 and @ we introduce another measure of 
localization. For each A^-dimensional eigenvector of the 
matrix Unm the information entropy is 

N 

J^(ui,...,ujv) =-^ Wnlnw„, ( 18 ) 

n—1 

where = |u„p, and = 1. We denote 



+ 1^ - ^ In + 0(l/iV), 

(19) 


where a = exp( 2 - 7 ) ~ 0-96, while ip is the digamma func¬ 
tion and 7 the Euler constant (~ 0.57721...). We thus 
define the entropy localization length Ih as 

Ih = N exp ■ (20) 

Indeed, for entirely extended eigenstates Ih = N. Thus, 
Ih can be calculated for every eigenstate individually. 
However, all eigenstates, while being quite different in 
detail, are exponentially localized, and thus statistically 
very similar. Therefore, in order to minimize the fluctu¬ 
ations one uses the mean localization length d = {Ih)} 
which is computed by averaging the entropy over all 
eigenvectors of the same matrix (or even over an ensem¬ 
ble of similar matrices of the same N but nearby k) 

d={lH)=Nexp{{J^H)-^P^°^)- ( 21 ) 

The localization parameter Pioc is then defined as 


^loc 


pL = 

N ~ N ' 


( 22 ) 


The parameter that determines the transition from weak 
to strong quantum chaos is neither the strength param¬ 
eter k nor the localization length £, but the ratio of the 
localization length £ to the size N of the system in mo¬ 
mentum p 



N\ ) 


(23) 


where C ~ tioc, the theoretical localization length 
Eq. (fT^ . was derived in Q. A is the scaling parameter 
of the system. The relationship of A to (3\oc is discussed 
in section VII of @ . 


III. THE DISTRIBUTION OF THE 
LOCALIZATION LENGTH AND OTHER 
LOCALIZATION MEASURES 


In this section we present the main results of the paper. 
First we demonstrate that the localization measures 2 /(t 


FIG. 1. We show {a) versus 2/{Ih) for matrices of dimension 
N = 3000, for 7 nearby values of k, namely k — ko ± jSk, 
where j = 0,1, 2, 3 and 5k = 0.00125, for kg = 3,4, 5,..., 19. 
The two empirical localization measures are clearly well de¬ 
fined, linearly related and thus equivalent. 


and Ih are very well defined, linearly related and thus 
equivalent. In Fig. [T]we show this in the diagram of the 
mean (cr) versus 2/{Ih), where both averagings are over 
all eigenfunctions for matrices of dimension N = 3000, 
for 7 nearby values of k around fco, namely k = ko ± 
jSk, where j = 0,1,2,3 and 6k = 0.00125, for ko = 
3,4,5,...,19. 

In the next Fig.[2]we show the relationship of the theo¬ 
retical C in Eq. (HU and the mean value of the empirical 
2/(cr) for ko = 3,4, 5,..., 19. It is clearly seen in Fig. [HJa) 
that there are strong fluctuations which we attribute to 
the fact that 2/a has a certain distribution with nonva¬ 
nishing variance, to be presented and described below, 
and that the theory of £ is too simple, as it corresponds 
only roughly to the value of 2/(cr). On the other hand, 
in Fig. (b) we see again that the two empirical local¬ 
ization measures are exactly linearly related. We should 
mention that in the cases of larger fc > 19 the slopes a are 
so small, and the localization too weak, that we cannot 
get reliable results, thus in this work we limit ourselves 
to the interval 3 < /c < 19. 

Thus we have demonstrated that the empirical local¬ 
ization measures are well defined, while the theoretical 
prediction for their mean values is not good enough. The 
reason is that the localization measures of a given fixed 
system (with fixed A = 10 and k) have a distribution 
with nonvanishing variance, which is out of the scope of 
current semiclassical theories, as they do not predict this 
distribution and the corresponding variance. This finding 
as the central result of the present paper is demonstrated 
in Fig. [21 The distributions are clearly seen to be close 
to a Gaussian, but cannot be exactly that as a is always 
a positive definite quantity. Its inverse, the localization 
length equal to 2/a, has a distribution whose empirical 
histograms are much further away from a Gaussian, so 
that in this sense a is the fundamental quantity. Indeed, 
as we will see, it corresponds to the finite time Lyapunov 
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FIG. 2. [Color online] (a) We show C versus ‘Ij{p) for matrices 
of dimension N = 1000 (crosses and solid fit line) and for 
matrices of dimension N = 3000 (stars and dashed fit line), for 
7 nearby values of fc, namely k = ko ± jdfe, where j = 0,1, 2, 3 
and 5k — 0.00125, for fco = 3,4, 5,..., 19. (b) We plot the 
mean value of 2/{N{a)) versus /3ioc for fco = 3, 4, 5,..., 19 
and 7 matrices of dimension N = 3000 with fc = fco ± jSk, 
where j = 0,1, 2, 3 and the step size Sk — 0.00125. 


exponent known in the theory of dynamical systems. 

As Ih and 2/(T are equivalent localization measures, 
the former one is expected also to have a distribution, 
which we demonstrate in the histograms of Fig. 21 

We have also analyzed how the localization measures 
vary in the semiclassical limit of the increasing value 
of the quantum parameter fc, at fixed classical param¬ 
eter K = 10. Indeed, the theoretical estimate of C in 
Eq. (1121) . at fixed K, and remembering k = K/t, shows 
that approximately the mean value of the localization 
length should increase quadratically with k, or equiva¬ 
lently, the slope a should decrease inversely quadratically 
with k. This prediction is observed, and is demonstrated 
in the Table HI and also in Fig. O It is also in agreement 
with the prediction based on the tight-binding approxi¬ 
mations in reference [Eq. (6)]. We give, in Table U 
the mean slope cr and the standard deviation of tr, as well 
as the mean value of the related quantity 2/1// and its 
standard deviation for various k = fco = 3,4, 5,..., 19, 
for each of them taking seven nearby values of fc, namely 
fc = fco ± jSk, where j = 0,1,2,3 and Sk = 0.00125, 
for matrices of dimension N = 3000. Each histogram 




o 



FIG. 3. We show the histograms of the slopes a for four sys¬ 
tems, matrices of dimension N = 3000, for each of them with 
seven different values of fc close to fco = 5, 9,13,17, namely 
fc = fco ± jhfe, where j = 0, 1,2,3 and 5k = 0.00125: (a) 
fco = 5, (b) fco = 9, (c) fco = 13 and (d) fco = 17. 
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FIG. 4. We show the histograms of Ih in (a) and 2/Ih in (b) for the system fc = 10 described by the matrices of dimension 
N = 3000. In both cases we show the Gaussian best fit. 



FIG. 5. We show log-log plots in (a) the mean slope (a) as a function of k, and in (b) the standard deviation of ct as a function 
of k. The htting by a straight line is only on the semiclassical interval 10 < A: < 19. In the former case the behavior is roughly 
as 1/fc^, in agreement with the theoretical estimate 1/k^ of Eq. (I12II . and in the latter case also like 1/k^, surely not as the 
theoretical estimate 1/k based on the Lyapunov exponents method in the reference [Eq. (9)]. 


for all ko was fitted with the Gaussian distribution and 
then the mean values and the standard deviations were 
extracted. All four quantities decrease to zero with in¬ 
creasing k, meaning that in the semiclassical limit the 
localization lengths monotonically increase to infinity, so 
that in this limit we have asymptotically extended states 
(no localization), and their standard deviation also goes 
to zero as 1/fc^, which is different from the tight-binding 
approximations in reference [Eq. (9)]. 

Next we want to study how does the distribution of 
the localization measure a behave as a function of the 
dimension N of the Izrailev model Eqs. (HMni). Since in 
the limit N ^ oo the model converges to the infinitely 
dimensional quantum kicked rotator, we would at first 
sight expect that following the Shepelyansky picture Q a 
should converge to its asymptotic value, which is sharply 
defined in the sense that the variance of the distribution 
of a goes to zero inversely with N. Namely, at fixed 
K and k Shepelyansky reduces the problem of calculat¬ 
ing the localization length to the problem of the finite 
time Lyapunov exponents of the approximate underly¬ 
ing finite dimensional Hamilton system with dimension 


2k. The localization length is then found to be equal to 
the inverse value of the smallest positive Lyapunov ex¬ 
ponent. In our case, the dimension of the matrices N of 
the Izrailev model plays the role of time. As it is known, 
and analyzed in detail in the next sections, the finite time 
Lyapunov exponents have a distribution, which is almost 
Gaussian, and its variance decays to zero inversely with 
time. Thus on the basis of this we would expect that the 
variance of a decays inversely with N. 

However, this is not what we observe. In the Table m 
we clearly see that at constant iL = 10 and A: = 10 the 
mean value of a is constant and obviously equal to its 
asymptotic value of = oo, while the variance of cr 
does not decrease with N, as 1/A, but is constant in¬ 
stead, independent of A. This is in disagreement with 
the banded-matrix models of the tight-binding approxi¬ 
mations and thus disagrees with the Eq. (9) of reference 
[i^l - and also disagrees with the Shepelyansky picture. 
The reason is that the associated Shepelyansky’s Hamil¬ 
ton system is only approximate construction, because 
with increasing A the matrix elements of the Floquet 
propagator (matrix) outside the diagonal band of width 











TABLE 1. 


The mean value and the standard deviation of the 
slopes a and 2jin as a function of 


N = 

A: = feo = 3, 4, 5,..., 19. For each k = ko we used 
= 7 X 3000 slopes a (see text). All quantities decay 
to zero in the semiclassical limit. 

A = 10 - 
k < a > 

N = 7 X 3000 (slopes) 
SDa < 2 /Ih > 

- N 

= 3000 {2/Ih) 

SD 2 /IH 

3 

0.06209 

0.01324 

0.062098 


0.01324 

4 

0.04327 

0.01073 

0.043272 


0.01073 

5 

0.04636 

0.00758 

0.046363 


0.00758 

6 

0.04030 

0.00974 

0.040303 


0.00974 

7 

0.04095 

0.00838 

0.040954 


0.00838 

8 

0.03004 

0.00756 

0.030047 


0.00756 

9 

0.03174 

0.00600 

0.031743 


0.00600 

10 

0.02835 

0.00539 

0.028355 


0.00539 

11 

0.02034 

0.00353 

0.020341 


0.00353 

12 

0.02014 

0.00321 

0.020143 


0.00321 

13 

0.01719 

0.0029 

0.017193 


0.00294 

14 

0.01750 

0.00289 

0.017509 


0.00289 

15 

0.01356 

0.00230 

0.013569 


0.00230 

16 

0.01221 

0.00194 

0.012213 


0.00194 

17 

0.00978 

0.00148 

0.009787 


0.00148 

18 

0.00855 

0.00128 

0.008550 


0.00128 

19 

0.00975 

0.00141 

0.009754 


0.00141 


TABLE 11. 

The mean value and the variance of the slope ct as a 
function of the matrix dimension N for a fixed 
system with AT = 10 and k = 10. Both are obviously 
constant. 


N 

A = 10 - fc = 10 

VaVa 

500 

0.102624 

0.00113224 

1000 

0.101170 

0.00112558 

2000 

0.100066 

0.00115575 

3000 

0.102217 

0.00110438 


2fc become important, and thus the dimension of the 
Hamilton system cannot be considered finite, constant 
and equal to 2fc, but increases with N. As a, consequence 
we have the constant value of the variance of a, and thus 
constant variance of the localization length C = 2/(t, and 
therefore the localization length has a distribution with 
nonvanishing variance even in the limit N = oo. This is 
precisely the reason why the semiclassical prediction of 
the localization length in Eq. m fails in detail and we 
find strong fluctuations in the plot of £ against the 2 ja of 
Fig. [5] The proper theory of the localization length must 
predict its distribution rather than just its approximate 
mean value. 


IV. NUMERICAL STUDY OF FINITE TIME 
LYAPUNOV EXPONENTS FOR THE 
CLASSICAL STANDARD MAP 

Finite time Lyapunov exponents of chaotic systems is 
a subject of not very much intense research. Taking an 
ensemble of uniformly distributed initial conditions of a 
uniformly chaotic (ergodic) system (with no islands of 
stability) we of course expect that for any finite time 
the Lyapunov exponents will have a certain distribution. 
With increasing time the mean value of each of them is 
expected to converge to the asymptotic Lyapunov expo¬ 
nent, and since the asymptotic Lyapunov exponent must 
be the same for all initial conditions, the distribution 
must converge to the Dirac delta distribution. Some early 
results on this topic go back to the 1980s, in the works 
of Fujisaka [i^, reviewed and summarized by Ott (4l| . 
Some details are not so important, as it turns out that the 
distribution becomes Gaussian very fast with increasing 
time, which we want to demonstrate in this section. 

In Fig. [S]we show the histograms of the positive finite 
time Lyapunov exponent for the standard map [Eq. ([3]) ] 
with K = 10, for the finite times (=number of itera¬ 
tions) t = 50,100,500,1000 in (a), (b), (c) and (d), re¬ 
spectively. The initial conditions, 200 x 200 on a grid, 
have been taken uniformly distributed over the square 
27r X 27r. Already at t = 50 the distribution is quite close 
to a Gaussian, and this trend increases very fast. At 
longer times like t = 2000,3000,4000,5000 it becomes a 
perfect Gaussian distribution (not shown). The variance 
decreases as 1/t, as it is demonstrated and analyzed in 

Sec. ED 


V. NUMERICAL STUDY OF FINITE TIME 
LYAPUNOV EXPONENTS FOR THE PRODUCT 
OF RANDOM SYMPLECTIC 2D MATRICES 

As it is well known the problem of quantum or dynam¬ 
ical localization is related to the Anderson localization 
model, within the framework of the tight-binding approx¬ 
imation, with hopping transitions between the nearest 
neighbors only. This goes back to the pioneering work of 
Fishman, Grempel and Prange [i^, as discussed in 0,0, 
and also reviewed in [d^. Assuming the nonvanishing 
nearest neighbor interaction only and the site disorder, 
the governing Schrbdinger equation is 0 

^n+l “f — Edn (^'^) 

where E is the eigenenergy of the eigenfunction, while 
is the fluctuating on-site potential, varying from site 
to site, with a certain probability distribution. Therefore 
we have the equation 
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MLE (standard map) 



FIG. 6. We show the histograms of the positive finite time 
Lyapunov exponents for the standard map [Eq. ([3])] with K = 
10 and times (number of iterations) t = 50,100, 500,1000 in 
(a), (b), (c) and (d) respectively. The initial conditions are 
on the grid 200 x 200 on the square [0, 27r) x [0, 27r). In all 
cases we show the Gaussian best fit. 


where the 2x2 transfer matrix T„ is given by 


T = 

TJ. - 


E-El -1 
1 0 


( 26 ) 


The determinant is equal to one, and W = E — El is 
drawn from a distribution, defined by a given model. 
Therefore the asymptotic properties of the eigenfunction 
coefficients a„ as a function of n are determined by the 
behavior of the product of the random transfer matri¬ 
ces, T = TnTn-i ■ ■ ■ T 2 T 1 . Everything is determined by 
the trace B = TrT. If \B\ > 2 the eigenvalues of T 
are real reciprocals, A > 1 and 1/A < 1. Typically A 
grows exponentially with n, and Mn = n~^ In A > 0 fluc¬ 
tuates with n, has certain distribution for each finite n, 
and the limit M = lim„_>.oo Mn exists. The latter is 
known as Furstenburg theorem [d^. Thus, for generic 
initial condition (ao,ai) the a„ will grow exponentially 
with n and only for a special initial condition, they will 
decrease exponentially with the rate Mn as n —>■ -l-oo, 
but still will increase exponentially in the backward di¬ 
rection n —^ — 00 . There are then exactly the eigenener- 
gies E for which the an decrease exponentially in both 
directions n —>■ ± 00 . In such case then Mn, the finite 
time Lyapunov exponent, is precisely the inverse value of 
the localization length in the n-space. Thus, for the finite 
system n < 00 , we shall have a certain distribution of the 
Lyapunov exponents M„. Indeed, this is observed in our 
numerical experiments shown in Fig. [71 for the box distri¬ 
bution of W, namely within the interval W G [—2,-|-2], 
for four values n = 50,100,500,1000, and for each of 
them for 10 000 realizations, drawn from the distribu¬ 
tion of W, and we see that the Gaussian approximation 
is very good, and becomes perfect for longer values of n, 
such as n = 2000,3000,4000, 5000 (not shown). 

We have also analyzed what happens if we replace the 
box distribution of W by other distributions, and con¬ 
vinced ourselves that the dependence on the details of 
the distribution of W is very weak, as the distribution 
of the finite time Lyapunov exponents is always Gaus¬ 
sian. In Fig. EKa) we show the histogram of the finite 
time Lyapunov exponents for n = 100 with the Gaussian 
distribution of W with zero mean and standard deviation 
equal to one. 

One might expect that things will be changed drasti¬ 
cally if the distribution of W is different, with diverging 
variance. In Fig.jHKb) we show the result for the Cauchy- 
Lorentz distribution of W defined as follows 

where b is the halfwidth at the half maximum, and we 
have chosen 6=1. We have taken the values inside 
the cut-off interval [—2, +2] and n = 100, and then the 
same thing for the interval [—100, -1-100] and n = 100 in 
Fig.EKc). We clearly see that the distribution is always 
Gaussian. 
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FIG. 7. We show the histograms of the positive finite time 
Lyapunov exponents for the product of random matrices 
[Eq. (E51l ] with W — E — E° uniformly distributed in a box 
W € [—2,+2], for n = 50,100,500,1000 in (a), (b), (c) and 
(d) respectively. In all cases we show the Gaussian best fit. 




MLE (random matrix — Cauchy w E [— 100 , 100 ]! 


FIG. 8. The histograms of the positive finite time Lyapunov 
exponents for the product of random matrices [Eq. (1261) 1 at 
n = 100 with (a) IF = E — E^ Gaussian distributed with 
zero mean and unit variance, (b) Cauchy-Lorentz distribution 
[Eq. (I27ll [ with W in the cut-off interval [—2, -|-2|, and (c) the 
same as (b) but W £ [-100,-1-100]. In all cases we show the 
Gaussian best fit which is excellent. 


VI. NUMERICAL STUDY OF THE DECAY OF 
THE VARIANCE OF THE DISTRIBUTION OF 
THE FINITE TIME LYAPUNOV EXPONENTS 

Finally, in this section we present numerical evidence 
for the theoretical expectation [dj that the finite time 
Lyapunov exponents have approximately Gaussian dis¬ 
tribution whose variance decreases inversely with time t 
(the number of iterations in the case of the standard map; 
Sec.HV]) and n in the case of the product of random ma- 
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FIG. 9. [Color online] The standard deviation of the positive 
finite time Lyapunov exponents for the standard map (stars) 
and for the product of random transfer matrices with the 
box distribution of W (empty boxes), as a function of time in 
log — log presentation, and their best fits. The slope is exactly 
- 1 / 2 . 

trices in the context of the unimodular transfer matrices 
of the tight-binding approximation to describe the An¬ 
derson localization, expounded in section El Indeed, the 
evidence is overwhelming, as shown in Fig. |9l where we 
plot the standard deviation as a function of time in log- 
log plot, showing that it decays inversely with the square 
root of time. 

In the context of our Izrailev model the dimension JV 
of the matrix plays the role of time. The width of the 
diagonal band is equal to 2k. Shepelyansky reduces the 
problem of the localization length to the determination 
of the smallest positive Lyapunov exponent (its inverse 
is the localization length) of the underlying finite dimen¬ 
sional Hamilton system with dimension 2k. Then, the 
finite time Lyapunov exponent should have some almost 
Gaussian distribution, whose mean tends to the asymp¬ 
totic Lyapunov exponent with iV —>■ oo and the variance 
should decrease to zero as 1/iV. 

If this picture were exact, then the mean localization 
length as a function of JV should converge to the asymp¬ 
totic value, which we do observe in Table IDof Sec. mil 
while the variance does not decay to zero, but rather 
remains constant, independent of N. From this we con¬ 
clude that even in the limit N ^ oo the localization 
length has a certain distribution with nonvanishing vari¬ 
ance, or more precisely, its inverse (the slope a) has an 
almost Gaussian distribution with nonvanishing variance. 
We believe that this is the cause of the strong fluctua¬ 
tions observed for example in Fig. Ha) of Sec. Uni The 
same observation applies to the scaling laws of the Figs. 9 
and 10 of our previous paper Q. 

VII. SUMMARY 

The main conclusion of this paper is the empirical fact 
based on our numerical computations of the eigenfunc¬ 


tions of the V-dimensional Izrailev model, that the local¬ 
ization length has a distribution with nonvanishing vari¬ 
ance not only for hnite N, but even in the limit N —>■ oo. 
This is the reason, we believe, for the strong fluctuations 
in the scaling laws which involve the empirical localiza¬ 
tion measures and the theoretical semiclassical value of 
the localization length. In the Shepelyansky picture Q 
this might seem to be a contradiction, but the resolu¬ 
tion of the puzzle is that in the limit of large N the 
finite dimensional Hamilton system extracted from the 
Floquet propagator of the quantum kicked rotator is not 
good enough, and therefore the matrix elements outside 
the main diagonal band of width 2k play a role, mak¬ 
ing the Hamilton system effectively infinite dimensional, 
with infinitely many Lyapunov exponents. This finding 
is a challenge for the improved semiclassical theory of 
the localization length, to derive and explain the discov¬ 
ered distribution function. On the other hand, the simple 
model of the Anderson localization based on the tight- 
binding approximation, with only the nearest neighbor 
interactions, described by the product of 2 x 2 unimodular 
matrices, has a finite dimension, as the transfer matrices 
are exactly two-dimensional, and therefore the variance 
vanishes in the limit of large dimensions as 1/n. The 
same conclusion applies to such a model with a hnite 
number of interactiM neighbors. Indeed, according to 
the references [i^, |4^ the variance of a should vanish 
as Var{a) oc l/(iVfc^), but our work shows that in the 
quantum kicked rotator this is not observed: the vari¬ 
ance does not depend on N, and decays with k faster 
than 1/fc^, namely as 1/fc^. Thus, here we found some 
important differences between the dynamical localization 
in the quantum kicked rotator and the Anderson tight- 
binding model of localization, and the Shepelyansky pic¬ 
ture, which rest upon the banded matrix models with 
hnite bandwidth. 

To summarize: We do not have yet a theory to de¬ 
scribe this behavior, namely the theory of the distribu¬ 
tion of the localization length, including the variance, 
rather than just its average value, as explained in the 
paper, but only the clear understanding of what is the 
reason for this behavior: The fact that the banded ma¬ 
trix model for the QKR is not good enough, one has to 
take into account also the (small but many) matrix ele¬ 
ments outside the main diagonal band, and therefore the 
Shepelyansky picture and approximation breaks down, 
meaning that the finite dimensional Hamiltonian system 
cannot capture the correct behaviour of the QKR. Thus, 
the problem is open for the future work. 
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